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Foil  bearings  are  intended  to  operate  at  high  speeds.  For  high 
speed  bearings,  in  general,  the  effects  of  fluid  inertia,  turbulent 
flow  and  heat  transfer  become  important  issues  in  analyzing  bearing 
performance  characteristics.  For  the  foil  bearings,  the  thermoelas¬ 
ticity  of  the  flexible  foils  adds  complications  to  the  analysis  of 
the  high  speed  bearings.  A  competent  advanced  analytical  tool  for  the 
foil  bearings  should  take  the  above-mentioned  issues  into  consideration. 

The  conventional  Reynolds  equation  approach,  for  several  decades,  has 
beer,  a  prevailing  method  for  hydrodynamic  lubrication  analysis.  The 
basic  assumptions  of  the  Reynolds  equation  a-re:  (1)  The  Reynolds  number 
is  ver>  small,  (2)  The  flow  across  the  film  is  insignificant.  And,  (3) 
the  fluid  is  isothermal.  Apparently,  these  assumptions  become  inadequate 
as  the  Reynolds  number  is  greater  than  the  order  of  magnitude  of  one. 

The^e  have  been  technical  papers  and  reports  contributed  to  the 
understanding  of  the  turbulent  lubrication,  the  inertial  and  the  thermal 
effects.  For  the  turbulent  fluid  film,  the  mi x i ng- 1 ength  [1]  and  the 
"wall  of  the  law"  [2]  and  [3]  are  the  two  most  popular  approaches. 
Launder  and  Leschriner  [4]  later  employed  the  low  Reynolds  number  k-e 
turbulence  model  [5]  to  treat  slider  bearings.  The  k-e  model  [6-7]  has 
beer,  proved  to  be  more  in  line  with  experimental  observations  than 
ether  turbulence  models  in  various  flow  situations  such  as  flat  plate 
boundary  layers,  pipe  flows  and  shear  flows  of  jets.  Although  the  U-c 
motel  adds  two  mere  equations,  the  t^i-bulence  kinetic  energy,  k,  and 
the  dissipation  rate  of  turbulent  energy,  e,  equations  to  the  computa¬ 
tion  process,  it  relates  the  turbulent  viscosity  to  the  lecal  turbulent 
energy  tn&nsport  and  presides  a  better  tool  tc  trace  the  history  cf  the 
turbulent  flow. 

In  dealing  with  the  fluid  inertia,  represented  by  the  convective  terms 
in  the  rro  ”  ertum  equations,  most  cf  previous  studies  integrated  the 
momentum  equations  across  the  film  gap,  a  process  similar  tc  the  deriva¬ 
tion  of  the  Reynolds  equation.  Since  the  velocity  prof i 'e  is  net  known 
beforehand  ir.  the  film  gap.  tc  carry  cut  the  '  r  teg'  at  icr ,  a  velocity 
profile  is  assumed  cr  Lu:  towed  •*  r  c  r.  the  nen-  1  ne  *"1  •  a"  solution. 


The  runner  surface  and  the  tear  ..j  ;  or  foil  )  surface  are  ret  parallel  . 
The  inclination  of  the  surfaces  results  in  a  transverse  velocity 
component  contributed  by  the  runner’s  tangential  velocity,  adding  to 
the  "low  across  the  film  gap.  This  transverse  velocity  component  becomes 
more  pronounced  as  the  runner's  speed  increases.  Consequently,  an 
assumed  velocity  profile  may  not  be  able  tc  reveal  the  local  flew 
behaviour  especially  for  the  high  speed  bearings.  The  acrcss-f i !m  inte¬ 
gration,  at  best,  can  only  catch  the  velocities  at  the-  two  ends. 

The  present  study  adepts  a  three  dimensional  Navier-Stokes  approach 
with  the  low  Reynolds  number  k-e  turb-'ence  mode!  by  Jones  and  Launder 
[3]  for  the  turbulent  flow  case.  The  Navier-Stokes  equations  are  written 
in  cylinderica!  coordinates  sc  that  the  curvature  effects  including  the 
centrifugal  and  the  coriclis  accelerations  can  be  taken  into  account. 

The  momentum  equations,  the  continuity  equation,  the  energy  equation, 
and  the  k-e  equations  are  solved  by  the  SIMPLE  ( Semi - Imp! i ci t-Pressuro- 
Li nked-Equat ions )  algorithm  developed  by  Patankar  and  Spalding  [8]  and 
elaborated  in  the  book  by  Patankar  [9]. 

The  thermohydrodynamic  solution  of  the  fluid  film  provides  the 
pressure  and  the  temperature  distributions  on  the  foil  surfaces.  The 
fluid  pressure  and  the  thermal  effect  cause  the  deformation  of  the 
foils.  Tc  analyze  the  thermoelasticity  of  the  foil,  a  thin  plate  model 
is  chosen  to  simulate  the  foil.  The  governing  thermoelasticity  equation, 
a  fourth  order  partial  differential  equation,  is  solved  by  finite- 
difference  approximations.  The  result  of  the  numerical  calculation  gives 
the  foil  deflection  which,  in  turn ,  changes  the  fluid  film  shape. 

The  interaction  between  the  the  rr.chy  dr  odynami  cs  of  the  fluid  film  and 
the  chermoelacti city  cf  the  foil  is  dealt  with  by  combining  the  two 
parts  of  computations  and  performing  iterations-  tc  reach  a  convergent 


Per  the  Phase  I  study,  two  computer  programs  are  developed,  the 
journal  bearing  and  the  thrust  bearing  programs.  The  two  kinds  cf 
bearings  have  quite  different  geometries.  Therefore,  the  equations 
for  the  two  bearings  are  not  the  same.  These  two  programs  can  be 
used  to  analyze  the  steady-state  performances  of  the  two  kinds  of 
bearings,  including  the  effects  of  the  inertia,  the  heat  transfer  and 
the  turbulent  flow.  The  programs  can  be  utilized  for  the  foil  bearing 
as  well  as  the  sclid-walled  bearing  analysis.  The  latter  is  simply  a 
special  case  as  far  as  the  program  capability  is  concerned. 

Computational  examples  are  presented  for  both  the  thrust  and  the 
journal  bearings.  The  performance  character i st i cs ,  such  as  the  load 
carrying  capacity,  the  frictional  coefficient  and  the  attitude  angle 
(  for  the  journal  bearing)  are  calculated.  Results  are  compared  with 
available  solid-walled  bearing  data,  [10]  and  [11]  to  verify  the 
accuracy  and  to  reveal  the  inertial,  the  thermal  and  the  turbulent  flow 
effects. 

The  Navier-Stokes  approach  differs  from  the  Reynolds  equation 
approach  not  only  in  solving  more  complicated  differential  equations 
and  realizing  the  above-mentioned  three  effects,  but  also  in  providing 
an  advanced  tool  for  visualizing  how  the  fluid  enter  the  bearing  and 
hew  the  flow  hehaves  locally.  The  continuation  of  this  study  such  as 
a  Phase  II  study  should  focus  on  the  trasient  case  and  the  combination 
cf  the  fluid  dynamics  and  the  reter  dynamics.  Hence,  a  better  picture 
of  the  dynamic  behaviors  cf  the  rote --bear  i  ng  can  be  obtained,  with 
limited  time  and  budget,  the  Phase  I  study  has  achieved  the  goal  cf 
establishing  the  fundaticn  for  pursuing  further  studies. 
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Basically,  the  fluid  film  in  either  a  journal  bearing  or  a  thrust 
bearing  is  equivalent  to  the  flow  between  two  inclined  surfaces.  For 
the  journal  bearing,  the  flow  is  between  two  eccentric  cylinders,  while 
for  the  thrust  bearing,  it  is  between  two  unparallel  disks.  To  describe 
the  three  dimensional  flew  behaviors  and  to  consider  the  curvature 
effects  including  the  centrifugal  and  the  coroil  is  forces,  cynlindrical 
coordinates  are  chosen  to  derive  the  governing  equations.  The  main 
geometrical  difference  between  the  journal  bearing  and  the  thrust 
bearing  lies  in  the  direction  of  the  film  gap;  for  the  former,  the  film 
gap  is  in  the  radial  (r)  direction,  while  for  the  latter,  it  is  in  the 
axial  (2)  direction.  Nevertheless,  both  bearings  have  a  common  feature 
that  is  a  main  flow  in  the  ci rcumferential  (0)  direction.  Due  to  this 
dominating  flow,  the  parabolization  concept  can  be  applied  to  the 
Nav 1 er-Stokes  equations  and  the  so-called  marching  technique  can  be 
employed  in  the  numerical  computation. 

Additionally,  because  of  the  eccentricity  of  the  cylinders,  in  the 
journal  bearing  case,  and  the  unparallel ness  of  the  disks,  in  the 
thrust  bearing  case,  the  main  flow,  V*  ,  has  components  in  the  other 
two  directions,  i.e.,  the  r  and  z  directions.  Furthermore,  the  shaft 
misalignment,  if  present  in  the  journal  bearing,  and  the  foil  deflec¬ 
tion  will  modify  fluxes  contributed  by  V$  to  the  r  and  z  directions. 
Consequently,  film-shape  fitted  coordinates  transformations  are  very 
much  desired  to  enable  the  realization  of  these  additional  fluxes 
in  both  the  r  and  z  directions. 

In  the  following,  formulations,  coordinate  transf ermat i ens  and 
numerical  solution  procedures  will  be  elaborated: 

A.  ~CRML»LATIQN£ 

The  Na\ i er-Stokes  equations  for  the  compress’ tie  flew  can  be  written 
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5>-^=f§.  -vp  +.mv\t  +  -jmV(V'U  )  (i 

The  continuity  equation  in  the  c>lindrical  coordinates  takes  the 
form 

If  4  i  TF  (?r\£>  4  T  ie  +-fe  )  =  0  (: 


The  energy  equation  in  terms  of  the  temperature,  T,  can  be  written 


fvU  +l/*^  =| £  +  U'VP  +  Vr  +  $  (2) 

where  R  is  the  thermal  conductivity  and  ^  is  the  dissipation  function 
which  can  be  expressed  as 

=  ( ^'*  V  )•£/  (*) 

where  <f>‘  -"S  the  stess  matrix  defined  as 

far  ^f\ 

<R.  VtJ  <5) 

Vir  Tl0  <7?J 

So  far,  the  abo\  e  equations  describe  the  th.ertr.oh>-  drotiynami  c  behaviors 
cf  the  fluid  film.  For  the  turbulent  flow,  ir,  the  high  speed  bearing 
case,  the  turbulence  disturbance  components  are  governed  by  the  Reynolds 
stress  equation: 

3  . r ,  ii..i  \  .  f  i  f  \rdUiU'3 

4ij«.  f  H  a,!,1 

4  -f  <  4  W)  -  (6) 
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where  i,j,k  denote  the  th’-ee  orthogonal  coordinate  directions,  and  U  is 
the  mean  flow  vector. 


It  is  prohibitive  to  solve  equation  (5),  even  numerically.  Several 
turbulence  models  have  been  developed  to  define  the  turbulent  viscosity. 
The  k-e  model  is  the  most  popular  and  well  tested  model  in  the  current 
stage  of  the  computational  fluid  dynamics. 

K  is  the  turbulent  kinetic  energy  and  e  is  the  rate  of  the  turbulent 
energy  dissipation.  They  are  defined  as 

k  =  "ir  Hi  U-i.  •  and  e  '  Vt  Tx%~  ^ 7 ' 

The  turbulent  viscosity  is  then  defined  as 

Qu?  K1-  /e  (8) 

Eased  on  Launder  and  Leschziner’s  simplification,  the  k-e  equations 
for  the  low  Reynolds  number  case  which  is  suitable  for  the  fluid  film 
in  lubrication  can  be  expressed  as 

■jO  KM*  *$£)  fjh  )  *MtK^-^  *  (f^)]  fe  -  -  0 

fr  tu<  +  If  j  +  ♦  (#?]  -  czf-^- 

♦  c3  <#f]  =  0 

Here,  it  should  be  noted  that  the  z-direction  is  the  direction 
across  the  film  gap;  for  this  case,  it  applys  tc  a  thrust  bearing. 

In  the  equations  (9)  and  (10),^  ,  ^,0"e  are  the  laminar  viscosity 
and  the  Prandtal  numbers  for  k  and  e.  C# ,  Cf,  C^and  C3  are  defined 
as  fell ows : 


(9) 


CIO) 
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Cl  =  1.44,  C2  =  1  .92  [  1  -  0.3  exp  (  -  Re*  )] 

C3  =  2.0,  <7k  =  1.C,  (7e  =  1.2 

Qu  =  0.09  exp  [  -  3 . 4/ (  1  +  Re  /50  )  ] 

Re*  =  f  k2  //x  e 

Re^  is  the  turbulent  Reynolds  number. 

Considering  the  steady  state  case  (  neglecting  -jfc-  terms),  applying 
the  parabolization  concept  (  neglecting  the  2nd  order  derivatives  with 
respect  to  0 ,  i . e . ,  -§5*  )  and  utilizing  the  continuity  equation,  all 
the  above-mentioned  governing  equations,  including  three  momentum 
equations,  two  turbulence  equations  (k-e)  and  one  energy  (temperature) 
equation,  can  be  written  in  a  generalized  form: 

-k  &  < rfVr<P )  +-\r  (PV/.*)  +  &  (Pt 


[convection  terms] 


[diffusion  terms]  [source  terms] 


(p  symbolizes  the  six  dependent  variables.  I^is  the  effective  viscosity 
Tf  =  M  +  represents  the  source  terms. 

The  convection  terms,  left  hand  side  of  the  equation  (11)  and  the 
diffusion  terms,  except!^,  are  standard  for  all  (p's.  The  Q  and 
for  all  s  are  tabulated  in  Table  1. 
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Tac'e  I  is  applicable  tz  the  th.-ust  beai'  ir.g .  F^r  the  journal  bearing, 
the  ro’es  of  r  and  z  are  exchanged.  Tc r.tc  cf  Table  I  can  be  modified 
azco:  di  ngl  >■  to  suit  the  journal  bear  *19  case. 

The  next  region  to  be  considered  is  the  foil  which  will  deform  under 
the  "rmlusnce  cf  the  fluid  film  pressure  and  the  thermal  effect.  The 
equat'icr.  describing  the  thermoelastic  behavior  cf  the  foil  which  is 
modelled  as  a  thin  plate,  is  adopted  from  Parl.us’  book  [12]: 

Vt  =  p  -  (  1  <*V*V  C2) 

In  writing  the  equation  (12),  the  deflection  caused  by  membrane 
forces  has  been  neglected.  This  is  based  or  the  assumption  that  the 
deflection  due  to  stretching  is  small  in  comparison  with  bending.  In  the 
equation  (12),  the  meanings  of  symbols  are: 


w  : 

deflection  in  the  direction  normal  to  the  fcil  surface 

E  : 

Young’s  modulus 

r  ■ 

Poisson’s  ratio 

t  : 

foil  thickness 

therma'  expans  :  or.  coef  *  .  s  i  ert 

P 

fluid  film  p j ecsui e 

*T: 

thermal  moment;  .  t  is  defined  ac 

rVi. 

1 

It 

J*.  \  T  ■*  d ■ 
t3  “ 

;  '  I 

The  factor  ]2a'-yi)  is  the  ^o-called  bending  st :  i -r  e  ss  .  Let 

• 

K 

E  t3 

rs 

ao-v1) 

For  bearings  with  combined  top  and  bottom  foils  (Fig.1)  such  as  MTI’s 
journal  and  thrust  bearings,  and  AiResearch's  thrust  bearing,  an 
equivalent  bending  stiffness  can  be  derived  as  follows: 


12  a -Ka) 

~s ~ (referred  to  [i2]  *  [i4]) 

and,  T  fa -  (14) 

t  U  -L.  if 


To  combine  the  equation  (12)  with  the  equation  (12),  we  take  a 
simplified  model  considering  one  dimensional  heat  conduction  across  the 
foil  thickness  with  boundary  conditions  depicted  in  Fig.  2. 


Fig.  2  Heat  Ccnduct’on  across  the  Foil 


Consequently,  the  temperature  distribution  across  the  foil  thickness, 
in  terms  of  the  wall  temperature,  Tw,  is: 


10 


Substituting  the  equation  (15)  into  the  equation  (13),  we  obtain 

2. 

yfl  li' t  [  •*  1  ^  A  A 

=  Tw  (Z  +  -~)  2  d2 


* Tw 


(16) 


Using  the  equation  (16),  the  thermoelasticity  equation  (12)  becomes 
K  Vw-  P  -  Tw  (17) 

For  the  journal  bearing,  p  and  Tw  are  functions  of  ©  and  z  while  for 
the  thrust  bearing  they  are  functions  of  r  and  0. 


B .  Coordinate  Transformations 

As  mentioned  previously,  due  to  the  film  thickness  variation,  the 
main  flow  contributes  fluxes  to  the  other  two  directions.  Since  the 
film  gap  direction  of  the  journal  bearing  differs  from  that  of  the 
thrust  bearing,  different  coordinate  transf ormat i ons  should  be  selected 
for  the  two  bearings. 

For  the  thrust  bearing,  we  choose  the  following  transformations  of 
coordinates  (referred  to  the  sketch  shown  in  Fig.  3): 


Fi 


y  •  w 


^reformations  (Thrust  Eeanng) 


Coordinate  Tr 


-Erik 

AR 

0  -  O.Vt 

& 

2  -Z,-n 

H 


(  18) 


g-fe-K) 

H 


3  -2c 

H 


where  A  R  =  rj.ut  -  r,rt  ,  the  difference  between  the  cuter  radius  and  the 
inner  radius.  6)  is  the  pad  angle  and  H  is  the  local  film  thickness. 

The  new  coordinates  represent  a  set  of  canonical  forms  which  define 
each  variable  in  a  unit  range ,  i.e.,  0  to  1 .  Applying  the  chain  rule  of 


f f erent 

iation,  we 
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dr 
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IF 
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dO 
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9  32_ 

aT  do 
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Substitution  of  the  coordinate  transformations  represented  by  the 
equation  (19)  into  the  general  transport  equation  (11)  results  in 

~y  :ijr  &  'rf’v'-<*’;i  +  t®“  w(  f 

*  #  vr+  $g- 

~r~  (IkJ  TF  (  *  "JJ*  !1"J*  3  *7 iff’  *sr 


(20) 


the  journal  bearing,  the  film 


■  e ,  r e f e r r . n g  to  r  1  g 
"c  ’  lowing  forms  : 


-  ,  we  select 


,  in  the  radial 
the  coord" rate 


c  i  rect i on . 

w  r*  £,  T.  S  "f  C.  T  T.  2;  “t.  *  C  t 


~i~  !-X-  -4- r  r94>  ( Vr  +  ilzHl  ^£L  w  x  (hr)  —L  )  ] } 

r  TT  ^FLrrVvVr  /  ae  ve  i  a*  T*;J' 

'  To"  W{^]  +  ~lT  JT[ 


~r  fjT  ~$F’(  +  iff 


]  +  it  ir{  %  i 


+  s. 


(22) 


'.r  the  equations  (20)  and  (23),  those  terms  representing  the  film 
cl. ness  variation,  namely,  ~r  and  in  (20),  and  ^|£*  and  i 


i  n 


(22),  confirm  the  fact  that  the  flow  in  the  c i rcumf erent i al  and  axial 
directions,  (or  ci rcumferential  and  radial  directions)  contribute  flux 
components  to  the  flew  across  the  film  gap.  For  high  speed  bearings,  the 
flux  component  contributed  by,  especially,  the  c i rcumferent i al  flew 
becomes  pronounced.  The  two  dimensional  approach  or  the  Reynolds  equation 
approach  may  net  be  able  to  reveal  locally  the  significance  of  this 
transverse  flux  increase. 


The  steady-state  continuity  equation,  for  the  thrust  bearing,  becomes 

~r~  ~5r  TT (  f  r  vr  '  *  "FiT  af"*  f  ^  ' 

i  ^rtf<V*+  ^  18-  If  VI!  =  0  :24! 

For  tie  journal  bearing,  the  continuity  equation  takes  the  form: 
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C .  hymeneal  Solution  Procedures 
; I )  Finite  Difference  Equations 

The  geometry  and  coordinates  of  the  journal  bearing  will  be  used  to 
illustrate  the  construction  of  finite  difference  equations. 


The  numerical  computation  will  be  carried  out  step  by  step,  the 
so-called  marching  process,  in  the  ci rcumferenti al  (©)  direction. 

T.e  range  of  ©,  such  as  the  arc  angle  of  a  partial  arc  journal  bearing, 
is  divided  into  a  finite  number  of  steps.  At  each  step,  the  r-z  plane 
is  discretized  into  a  grid  network.  Fig.  5  shows  a  central  point  P 
surrounded  bye  four  points  denoted  as  W  (west),  E  (east),  5  (south) 
and  N  (north).  Then,  a  control  volume  is  constructed  by  selecting  the 
ndd'e  grid  points,  the  lower  case  w,e,s  and  n,  on  the  <--z  plane,  ar.d 
a  dept:’:  49.  The  shaded  area  in  Fig.  5  indicates  one  of  the  control 
vc' ume  surfaces  on  the  r-z  plane. 

Tie  general  transport  equation  and  the  continuity  equation  a^e 
ntegrated  over  the  control  volume  to  obtain  the  finite-difference 
equations.  ~he  transport  equation  represents  a  convection-diffusion 


Ir.  computat i cnal  fluid  dynamics  (CFD),  it  is  well  known  that  the 
commonly  used  central -d i fference  scheme  will  lead  to  unstable  solutions 
as  the  convection  is  much  stronger  than  the  diffusion.  To  prevent 
numerical  instability,  the  upwind  difference  scheme  will  be  used. 

To  simplify  the  expression  of  the  transport  equation,  we  define 
the  following  symbols: 

F,  =  f  V* 
f2  =  f  ve 

F*  =  ?Vr+  *6  +  H-r)  |f  fV2 

/on 

D|  -  r^c, 

=  {i  +  [(i-r)-ff  hr+ 

Ar  =  Lr®  (92-  ©,  )  (2fi-  %,) 

A  a  =  LH  ( -  hj  )  ( Zg  -  Z*r) 

i  1  aj  a2  r  z' 

Al  =  *  H  (r"  "  rS  }  (eZ~  G*  } 

Furthermore,  subscripts  n,s,e,w  and  up  are  used  to  indicate  five 
locations;  up  is  the  up-stream  plane,  i.e.,  £t  plane. 

The  upwind  difference  scheme  can  be  expressed  as 


<  F.  <t>  )e  -  5  r,e  •  c  8  “  ^  B  -  r,e  ,  C  5 


here  3 

a,  t  fi  =  maximum 

^  ^  '  c,  K  ' 

w  .  <  a  ,  O  > 

-  C  c 

ie 

■;  0,  then  (=>;e 

=  p.e<£P 

*  *  le 

*.  Then  :  r,  <p  ;e 

=  -«e  ^ 

c 

he  abcve-def i ned  symbols,  the  transport  equation  becomes 


ap^>=  ae4fc+  a*<k  +  +  as^s*  £c/*  +  Fiup  ^pA©  :28' 

where 

a£  =  (5  -F.e.  0  5  +  rrfr}  Aae 

£iV  =  (0  F.«,  0  |  +  2^-}  Aj* 

a*  =  {0  *FJrt.  0  I  +  -jfeft-J  Ar„  (29) 

as  =  (0  FiS,  o  8  +  Ars 

ap  =  a£  +  aW  +  atJ  +  as  +  ( F‘«  -  fiur  )  +  (f>n  -  F,S  )  +  Sp4 
The  source  term  has  been  split  into  two  parts: 


The  continuity  equation  becomes 

'^trt  Ar#»  "  *>s  Ars  )  +  (F»e  Aae  "  Fu*r  ) 

-  (Fj,  dr.  A01-  Fa,  up  A®,  )  =  0  (31) 

c  -ite  difference  equations  for  the  the  mice  1 act i c i ty  equation  (17) 
arc  to  be  constructed  on  the  r-6  plaice  for  the  thrust  bearing  and  on 
the  z-r  plane  for  the  journal  tearing.  Here,  we  take  the  journal 
oear-.rg  case  as  an  illustration  example.  The  same  derivaticr  procedures 
3’"t  applicable  to  the  thrust  bearing  case. 


We  set  up  grid  networl.  cr  the  bearing  foil  and  select  eight  grid 
points  surrounding  the  central  point  P.  The  governing  equation  has 
bi-harmcnic  functions;  it  needs  9  points  to  approximate  the  fourth 
derivatives  with  the  use  of  central-differencing  scheme  whi-ch  is 
suitabel  for  this  case.  '  i 

Referring  to  the  grid  system  shown  in  Pig.  5,  the  thermcel ast i c i ty 
equation  becomes 
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f  ~  2  •  5  2nd  System  f:r  Thermcel  ast  i  c  i  ty  Eq. 
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dures  of  SIMPLE  Algorithm 


D.  Frcce 

"he  transport  and  continuity  equations  are  to  be  solved  by  the 
SIMPLE  algorithm  described  in  Patar.kar’s  bock  [9].  The  word  SIMPLE 
stands  for  Semi-Implicit  Method  for  Pressure-Li nked  Equations. 

The  major  operations  of  the  SIMPLE  algorithm  can  be  summarized  as 

' 1 )  Guess  the  pressure  field  denoted  as  p*. 

(2)  Solve  the  momentum  equations  to  obtain  V*  ,  v£  and  Vz. 

(2)  Solve  the  pressure-correction  equation  to  obtain  p’. 

(4)  Revise  the  pressure  p  by  adding  p’  to  p*. 

(5)  Correct  the  velocity  components  by  using  p’. 

(6)  Solve  the  transport  eqution  for  other  ^’s  such  as  temperature,  T 
k  and  e. 

Referring  to  the  finite  difference  equation  (28),  the  momentum 
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Let 


p  =  p*  +  p’ 

V©  =  Ve  +  vi 

•  (24) 

Vf  —  Vf  +  Vf 

\fe  =  V*  +  v* 

The  primed  terms  p’,  Ve’  ,  Vr’  and  Va  are  correction  terms. 

Substitution  of  the  equation  (34)  into  the  equation  (33)  results  in 

ve  =  s XmJj  +  (  Cj  -  Pd'  )  0®S(P<;  -  P;  )  o® 

v)  =  S*ivri  ♦  <  R,'  -  Pp  )  -  Pj,'  )  Dp  (35) 

vi  =  +  (  PL  -  «v'  )  D?~<K;  -  Dp  !  D* 

The  drop  of  the  terms  of  accounts  for  the  name  of 

Semi-Implicit  method.  Since  the  equation  (35)  has  to  satisfy  the 
continuity  equation  and  the  above-mentioned  operations  are  to  be 
iterated  until  a  converged  solution  is  reached,  the  Semi-Implicit 
process  will  not  sacrifice  the  accuracy  of  the  calculation. 


Substitution  of  the  equation  (25)  in 
results  in  a  pressure  correction  equat 

“P  pp  ”  aE  KE  ~‘iv  PW  aN 


+ 


o  the  continuity 
on  which  has  the 

aS  ps  +  b  +  au 


equation  (31) 
following  form: 
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and , 

b  =  ft,  Arn  (Vr;  +  B2m  V<£  +  Bm  Van  J  -ft  ( V*  +  Bas  V*s  +  Bls  v£  ) 

+  fd  vet)  Agj  "  iuveu  Aca  +je'4<  A*<  “  furvzur  A  aw 

With  the  establishment  of  the  quations  (34)  and  (26),  the  SIMPLE 
algorithm  procedures  can  be  readily  carried  out.  The  only  step  left  ou 
is  the  iteration  process.  For  the  present  case,  either  the  journal 
bearing  or  the  thrust  bearing,  the  dominant  flow  is  in  the  circum¬ 
ferential  (©)  direction;  especially,  in  the  journal  bearing,  the  flew 
is  in  cyclical  pattern.  Consequently,  the  numerical  calculation  is 
marched  in  the  9  direction  and  iterated  in‘0  until  a  converged 
solution  is  obtained. 

The  solutions  of  the  transport  equation  provide  the  pressure  and  the 
temperature  distributions  on  the  foil  surface.  Therefore,  the  numerics 
calculation  of  the  thermcel astici ty  equation  (32)  can  be  solved  by  any 
a'gebraic  equation  solver.  The  Gauss-Seidel  (with  over-relaxation) 
iteration  is  employed  to  carry  cut  the  calculation  of  the  foil 
deflection. 


*'•  c  foil  deflection  changes  the  fluid  film  shape,  which  in  turn 
''cuts  the  momentum  and  the  temperature  fields  in  the  fluid  film,, 
h  i  ■:  '  1  u d-"oi '  interaction  has  tc  be  iterated  to  reach  a  converged 
tatus  with  a  constraint  such  as  that  the  supporting  lead  ic  aqua1  • 
h„-.~  t  .ve  ght  or  the  fr  ct-cna7  f  c-'-.e  is  equal  tc  •  squired  .a' us. 


RESULTS  AND  DISCUSSION 


Load  capacity  (dimensional  and  dimensionless),  friction  force 
(dimensional  and  dimensionless)  and  attitude  angle  (for  the  journal 
bearing)  are  calculated  as  follows: 


(I)  Journal  Bearing 

Wi  =  W  sin^  =  f  (pw  -  p&  )(sin@)  rs  d0  dz 

W2.  =  W  cos^  =  ft  fZJpw  -  px)(cos©)  rs  d©  dz 
Jo  •'O 

where  W  is  the  dimensional  load  capacity; 

p^,  is  the  pressure  distribution  on  the  shaft  surface; 

is  the  ambient  pressure; 

£  is  the  length  of  the  shaft. 


Hence,  _________ 

W  =  /w i  +  W l 


The  attitude  angle  is 


The  dimensionless  load  capacity  is  defined  (to  be  compatible  with 
the  definition  given  in  [11])  as: 


W 


W 

D 


The  friction  forces  are  : 

f 

Fs  3  JJ.%,ias  aa 

f»  =  t  r‘  d® az 


(on  the  shaft  surface) 
(on  the  bearing  surface) 


where  is  the  shear  stress  on  the  shaft  surface- whi  le  Z^Ais  the 


shear  stress  on  the  bearing  (foil)  surface. 

The  dimensionless  friction  coefficient  is  defined  (referred  tc 
as 


[11]) 


.7  = 


px4dw 
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(II)  Thrust  Rearing 


The  load  capacity  is  calculated  by  integrating  the  pressure  over 
the  area  of  the  pad. 


Wl  =  -C£(fk"  V  r  d®  dz 

The  total  load  capacity  is  the  summation  of  the  loads  on  all  pads. 
W  =  N(pad)*Wi 

Dimensionless  load  capacity  is  defined  (referred  to  [10])  as: 

W  =  -fi-  JCL- 

*  PaLB 

where  B  =  Rt-  R,  and  L  =  /3  *(Ra  +  Ri)/2. 

The  friction  force  is  calculated  by  integrating  the  shear  stress  ; 
the  bearing  surface  over  the  area  of  the  pad. 


-  •£  fa' 


dr  d9 


The  dimensionless  friction  coefficient  is  defined  as: 


\ALHx 

where  -A  is  the  compressbi 1 i ty  number  and  Ha  is  the  minimum  film  ga 

Computations  of  numerical  examples  are  carried  out.  The  results  c 
the  computations  are  presented  in  the  following: 

(A)  Journal  Bearing 

"or  the  journal  bearing,  we  consider  a  bearing  and  a  shafts  with, 
input  data: 

RS(shaft  radius)=  0.0254  meter 
RB(bearing  radius)=  0.02542673  m 
RPMs  5CCOO.O 

ALEN(bearing  1enght)=  0.0508  m 
PAMB(ambient  pressure):  101353.0  newton/sq.  m 
HC(radial  clearance)=  0.000026726  m 
E£(radial  displacement):  0.000015036  m 


( 


BETA( arc  angle)=  6.2832  rad.  (full  journal  bearing) 

ALMULAM( 1  ami nar  viscosity)=  0.00002143  kg/m-s  (air) 

DEN(density )=  0.9611  kg/m 

Hence,  > 

_A.  (compressibility  number)=  6.0 
L/D( length  to  diameter  ratio)=  1.0 
&  (eccentricity  ratio)=  0.6 

The  computed  results  in  terms  of  dimensionless  load  capacity,  friction 
coefficient  and  attitude  angle  are  compared  with  the  data  given  by 
Wu  [11]  using  the  Reynolds  equation  approach,  namely,  isothermal  flew 
without  the  effects  of  inertia  and  turbulence. 


3CN8BRGJ 

Reynolds  Eqr, .  Approach 

Isothermal 

Non- 

i sothermal 

w 

1.137 

1  .157 

1 .023 

<$> 

19.5 

r 

- 

3.93 

(or.  bearing 
3.45 

(on  shaft) 

5.96 

(on  bearing) 
4.68 

(on  shaft) 

3.9 

t 

The  comparison  is  made  for  laminar  flow.  The  load  capacity  obtained 

by  the  Navier-Stokes  approach  is  higher  than  that  obtained  by  the  » 

Reynolds  equation  solution.  The  non- i sothermal  case  carries  higher 

load  than  the  isothermal  case.  The  fluid  inertia  enhances  the  load 

capacity,  and  also  increases  the  friction  slightly.  On  the  other  hand, 

the  temperature  rise  causes  the  increase  of  viscosity  Of  gas,  and 

contributes  to  the  increase  of  load  capacity.  However,  the  temperature 

change  also  increases  the  friction  considerably . 
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To  simulate  a  foil  journal  bearing,  we  consider  a  full  journal  bearing 
with  a  singTe  foil  leaf,  combination  of  the  top  and  the  bottom  foils, 
which  is  similar  to  a  MTI’s  bearing.  The  foil  deforms  under  the 
influence  of  the  hydrodynamic  pressure  and  the  heat  transfer  resulted 
from  temperature  variation.  The  foil  is  given  an  allowable  maximum  sway 
1  displacement,  SWYMX=  0.00005  meter;  the  rest  of  data  remain  the  same 

as  the  previous  example.  The  variations  of  load  capacity  and  friction 
*  coefficient  are  then  observed. 

The  comparison  of  the  solid-walled  bearing  with  the  foil  tearing 
is  tabulated  in  the  followint: 


Solid-Walled  Bearing 

Foil  Bearing 

Isothermal  1.137 

4 

.  «  w  W  w 

w 

Non-  1.157 

isothermal 

1  .323 

4> 

Isothermal  12.8 

12.0 

Non-  13.2 

i sothermal 

i  *  n 

V 

Isothermal  3.33 

(on  bearing) 

3.45 

(on  shaft) 

o  A*y 

•  *♦  w 

(or,  bearing) 
2.43 

(on  shaft) 

Non-  5.96 

isothermal  (on  bearing) 

4.68 

(on  shaft) 

7.48 

(on  bearing) 
n  o 

W  •  —  W 

(on  shaft) 

The  deformation  of  the  foil  cause  the  slight  increase  of  load 
capacity.  The  drop  of  friction  in  the  isothermal  case  is  due  to  the 
widening  of  the  film  gap.  For  the  non-isotherma'  case,  the  local  thermo¬ 
elastic  deformation  of  the  foil  may  create  pockets.  Consequently,  the 
friction  coefficient  changes  considerably. 
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For  the  case  of  the  turbulent  flow  in  the  journal  bearing,  the  data 
of  the  previous  example  are  modified  to  have: 


RPM=  52000 
RS=  0.C508  m 
RB=  0.0509326  m 
hC=  0.00013363  m 
ES=  0.00008018  m 
ALEN  =  C. 1016  m 


The  L/D  ratio  and  the  eccentricity  ration  remain  the  same.  The  Reynolds 
number  increases  to  be  1658.  The  bearing  compressibility  number  becomes 
1.  Results  comparing  with  the  Reynolds  solution  are  tabulated  in  the 
fol lowing : 


\ 


3DNSBRGJ 
sol i d  wal 1 

3DNSBRGJ 

foil 

REYNOLDS  SOL. 
solid  wall 

w 

0.409 

0.490 

0.20 

n 

4.98 

(on  shaft) 

3.31 

(on  shaft) 

2.5 

(on  shaft) 

The  decrease  of  friction  for  the  foil  bearing  may  due  to  the  increase 
of  the  film  gap. 

(B)  Thrust  Bearing 

We  consider  a  thrust  bearing  with  the  following  input  data: 

RPM=  3500C 

PAM8=  "01253  N(newton )/sq . m 

P.1:  0.1  in  r 

R2=  0.2  m 

BETA=  0.667  rad.  * 

H 1 =  0.00016703  m 
H2=  0.000012527  m 

AMULAN  (viscosity)=  2.143E-5  kg/m-sec 
DEN( densi ty ) =0 . 96 1  kg/cube  m 
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as : 


The  bearing  compressibility  number-/!,  L/B  and  H1/H2  are  calculated 


The  computed  results  for  the  isothermal  case  are  compared  with  the 
Reynolds  equation  solution  for  a  sector  thrust  bearing  giver,  in  [15] 
and  [10].  The  results  in  terms  of  load  capacity  and  friction  coefficient 
are  tabulated  in  the  following: 


The  dimensionless  load  capacity  by  the  Navier-Stokes  solution  is 
higher  than  that  by  the  Reynolds  equation  solution.  The  reason  is  that 
the  bearing  speed  is  high  and  the  inertial  acceleration  is  si 
The  friction  force  is  ottained  by  selecting  3-pcint  velocity 
near  the  wall  for  the  shear  stress  and  integrating  over  the  bearing 
surface.  The  great  difference  of  friction  coefficient  may  be  partia'ly 
due  to  the  coarse  grid  in  the  numerical  computation.  However,  the 
increase  of  friction  force  by  the  N-S  approach  is  physically  sc-nd. 
With  the  deflection  of  the  foil,  the  load  capacity  increases  and  the 
friction  coefficient  decreases  slightly. 
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PROGRAM 


(  Two  computer  codes,  3DNSBRGJ  and  3DNSBRGT,  are  developed.  The  names 

stand  for  three  (3)  Dimensional  Navier-Stokes  Bearing,  Journal  bearing 
„  and  Thrust  bearing  programs.  As  discussed  previously,  the  geometrical 
configurations  of  the  two  bearings,  the  journal  bearing  and  the  thrust 
bearing,  are  different.  The  coordinate  transf ormati ons  are  different. 
The  fina’  forms  of  the  equations  are  not  the  same  for  the  two  bearings 
Hence,  it  is  more  efficient  to  perform  the  calculations  in  separate 
codes  than  in  a  single  general  code.  Nevertheless,  both  codes  have  the 
same  structure  which  consists  of  subroutines  performing  designated 
functions : 

Subroutine  Function 


DATAIN 
INITIAL 
UPSTREM 
VISCOS  1  GAMA 

FILM7K 

FLUX 

BOUNCRY 

SOURCE 

COEFFSL 

EWPNS 

COEFMOD 

SOLVE 


1  WLOAD 

ELASTHM 

*  CHECK  &  CHECKA 


accepts  input  data 
defines  initial  conditions 
establishes  upstream  conditions 
calculate  the  effective  viscosity 
coefficients 

determines  film  thickness 
calculates  fluxes 
set  up  boundary  conditions 
calculates  source  terms 
calculates  transport  and  pressure 
correction  equations 

calculates  coefficients  of  a'gebraic  eqr.s. 
modifies  coefficients  in  algebraic  eqr.s. 
solver  for  algebraic  eqns.,  using  both 
block  tri-diagonal  matrix  and  Gauss-Seidel 
plus  over-rel axatior.  interaticr. 
calculates  performance  charateri sti cs 
performs  thermoel ast i ci ty  calculations 
print  out  value  of  variable  arrays  such  as 
density,  viscosity,  velocity  components, 
pressure  and  temperature. 


29 


The  MAIN  program  controls  the  step-by-step  marching  process.  It  also 
writes  the  results  of  the  current  iteration  on  a  disk  file  ar.d  use  this 
disk  data  as  the  initial  data  for  the  next  iteration.  In  addition,  it 
writes  data  on  disk  files  for  later  calculations  of  bearing  performnace 
charateristics  such  as  load  capacity,  frictional  force  and  dimensionless 
parameters. 

The  program  reads  input  data  in  the  format  of  NAMELIST.  In  the 
following  all  NAMELISTs  for  the  journal  bearing  and  the  thrust  bearing 
will  be  explained  respectively: 

(I)  Journal  Bearing 


NAMELIST 

VARIABLE 

MEANING 

3RGDAT 

RPM 

rpm  of  shaft 

PAMB 

ambient  pressure 

RS 

radius  of  shaft 

RB 

radius  of  bearing  bore 

BETA 

subtended  angle  of  foil 

BETAS 

starting  arc  angle  of  foil 

BETA1 

ending  arc  angle  of  fluid  film 
on  foil 

HC 

radial  clearance 

EE 

eccentric  displacement  of  shaf 

ALEN 

length  of  shaft 

FLUIDAT 

AMU LAM 

laminar  flow  viscosity 

PR  LAM 

laminar  flow  Prandtl  number 

DEN 

initial  fluid  density 

CP 

heat  capacity  at  constant 
prssure 

CPDCV 

ratio  of  Cp  to  Cv 

GASCON 

gas  constant  of  fluid 

PR 

Prandtl  numbers  for  al 1 
dependent  variables 

30 


AKFAC 


WM 

(  PIN 

TIN 

PCRIT 

* 

GRIDAT  A1,A2,A3 

B1 ,B2,B3 

N 1  ,  N2 , N3 
Ml , M2 , M3 

VOLDAT  L,M 

IMAX , JMAX 

PFAC 

IPF 

PFA.PFB 

FMEG 

CCNTL  NSWP 

IXY.ISWP, JSWP 


turbulent  flow  constant  for 
initial  condition  (4.3E-6) 
molecular  weight  of  fluid 
initial  pressure 
initial  temperature 
pressure  correction  equation 
criterion  (2.0) 

3  divisions  of  computation 

domain  in  ©  direction 

e.g.  0.25,0.5,0.25 

same  as  A1,A2,A3  but  in  z 

direction,  e.g.  0.2333,0.3323, 

0.3333 

number  of  grid  points  in  A1,A2,A3 

divisions  respectively 

number  of  grid  points  in  31,22,23 

divisions  respectively 

number  of  control  volumes  in  :  i 

r  directions 

total  grid  points  in  z  and  r 
directions;  IMAX=L+1 ,  JMAX=M+1 
relaxation  factor  on  pressure 
correction  (0.95) 
initial  pressure  profile 
selector  0  for  sin  &  cos  profi'e 
1  for  sin  profile  only 
initial  pressure  profile 
ampl i tude 

over-relaxation  factor  for 
all  dependen  variables( 1 . 5-1 .£) 
tri-diagonal  block  matrix  solver 
sweep  number  for  each  variable 
(8*2) 

tri-diagonal  block  matrix  sweep 
index  (1,1,1) 
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ISOLV 


ELSCON 


NNV 


NITER 


(II)  Thrust  Bearing 


BRGOAT 


on  and  off  key  for  solving  the 
transport  eqn.  for  each  dependent 
variable;  1  solve  the  eqn.  C  by 

pass 

dependent  variables  are  in  the 
order  of  Vr(1 )  ,Vz(2)  ,Ve  (3)  ,p’  (4) 
HH(total  enthalpy)(5) ,  k(5),e(7) 
T(8) 

total  number  of  dependent 

variable,  NNV=8 

total  number  of  momentum  eqn. 

iteration 


NELS 

total  number  of  thermoelasticity 
eqn.  iteration 

REYCRT 

critical  Reynolds  number 

ALPHA 

thermal  expansion  coefficient 
of  foil  material 

PATHON 

Poisson’s  ratio  of  foil  material 

YN 

Young’s  modulus  of  foil  material 

T1 

top  foil  thickness 

T2 

bottom  foil  thickness 

WAVL 

half  wave  length  of  bump (bottom) 
foil 

SWYMX 

maximun  displacement  of  foil  sway 

RPM 

rpm  of  shaft 

PAMB 

ambient  pressure 

R1 

inner  radius  of  disk 

R2 

outer  radius  of  disk 

BETA 

subtended  pad  angle 

BETA1 

subtended  fluid  film  angle 

HI 

maximum  film  gap 

H2 

minimum  film  gap 

HMEAN 

mean  film  gap 
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FLUIDAT 


AMU LAM 
PRIAM 
DEN 
CP 

CPDCV 

GASCON 
PR 

AKFAC 
WM 
PIN 
TIN 
PCRIT 

GRIDAT  A1.A2.A3 

B1 , B2 , B3 

N1 ,N2 ,N3 
Ml , M2 , M3 

VOLDAT  L.M 

IMAX, JMAX 

PFAC 

IPF 

PFA.PF3 

CONTL  NSWP 

IXY.ISWP, JSWP 
ISOLV ( 8 ) 

NITER , NELS 
REYCRT 

ELSCON  ALPHA 

PATHON , YN 
T1  ,T2 


laminar  flow  viscosity 
laminar  flow  Prandtl  number 
initial  density  of  fluid 
heat  capacity  at  constant 
pressure 

heat  capacity  ratio  of 
Cp  to  Cv 

gas  constant  of  fluid 
array  of  Prandtl  numbers  fc 
all  dependent  variables 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing  but 
in  r-di recti  on 
same  as  journal  bearing 
same  as  journal  bearing 
number  of  control  volumes  i 
r  and  z  directions 
total  grid  points  in  r  and 
di rections 

same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  bearing 
same  as  journal  .bearing 
same  as  journal  bearing 
top  and  bottom  foil 
thicknesses 


Typical  input  files  for  both  the  journal  and  the  thrust  bearings  are 
provided  in  the  following: 

(I)  Journal  Bearing  Input  Data  File: 

$BRGDAT  RPM=52000 . 0 , PAMB=1 0 1352.97,RS=5.08E-2,RB=5.093363E-2, BETA=6 . 2322 ,  ' 

BETAS=  0.0, BETA1 =6 . 2832 ,HC= 1 .3363E-4 , EE=8 . 01 73E-5 , ALEN=1 .016E-1 ,  $END 
$FLUIDAT  AMULAM=2. 1 43E-5 , PRLAM=0 . 7 , DEN=0 . 96 1 1 , CP= 1 009 . 0 , CPDCV= 1 .4,  * 

GASC0N=8 . 3 1 4E+3 , PR ( 1 )=1. 0.1. 0,1 . 0,1. 0,0. 9, 1.0,1 . 3 , 0 . 9 , 0 . 9 , AKFAC=4 . 2E-7 , 
WM=23.93,PIN=101352.97,TIN=367,PCRIT=2.0,  SEND 
$GRIDAT  A1=Q. 25, A2 =0.5, A3 =0. 25, B1=0. 3333, B2 =0.3333, 33= 0. 3333, N1=7, 
N2=10,N3=7,M1 =4 ,M2=5 ,M3=4 ,  $END 

$VOLDAT  L= 1 0 , M= 1 0 , IMAX= 1 1 , JMAX= 1 1 , PFAC=0 . 75 , IPF=0 , PFA=0 . 55 , PFB=0 .55, 

FMEG= 1.0,  $END 

$C0NTL  NSWP( 1 )=2 , 7*2 , IXY= 1 ,ISWP=1 , JSWP=1 ,IS0LV( 1 )=4*1 , 3*0 ,0 ,NNV=8 , 

NITER= 1 ,NELS=2,REYCRT= 1000.0,  $END 
SELSCON  ALPHA=3 . 3333E-6 , PATHON=0 .4091 , YN=2 . 1374E+1 1 ,T1 =0 . 0001 524 , 

T2=0. 0001016 , WAVL=0 . 005 , SWYMX=5 . OE-5 ,  SEND 

(II)  Thrust  Bearing  Input  Data  File: 


SBRGDAT  RPM=35000 . 0 , PAMB= 101352. 97, R1=0. 1 , R2=0 . 2 , BETA=0 . 667 , 
BETA1 =  0. 667, H1=0. 00016703, H2=8. 3515E-5, HMEAN=1 .2527E-4,  SEND 


SFLUIDAT  AMULAM=2. 143E-5 
GASCON  =  3 . 3 1 4E+3 , PR( 1  )  =  1 


, PRLAM=0 . 7 , DEN=0 .9611 , CP=1 009 . 0 , CPDCV= 1 .4, 
.0,1. 0,1. 0,1. 0,0. 9, 1.0, 1.3, 0.3, 0.9, AKFAC=4 


WM=28 . 93 , PIN=1 01 352 . 97 ,TIN=367 , PCRIT=2 . 0 ,  SEND 
SGRIDAT  A1=0. 25, A2=0. 5, A3 =0.25,81=0. 3333, B2=0. 3233, 33=0. 3332, N 1=7, 
N2=10,N3=7,M1=4,M2=5,M3=4,  SEND 


SVOLDAT  L= 1 0 , M=6 , IMAX=1 1 , JMAX=7 , PFAC=0 . 99 , IPF=0 , PFA=C . 25 , PFB=0 . 35 , 
SCONTL  NSWP( 1 ) =2 , 7*2 , IXY= 1 , ISWP=1 , JSWP=1 ,ISCLV( 1 )=4*1 , 3*0,0,NMV=S, 
NITER= 1 , NELS=2 , REYCRT=5000 . 0 ,  SEND 
SELSCON  ALPHA=3 . 333E-6 , PATHON  =  0 . 409 1 ,SWYMX=1 . E-5 , YN=2 . 1374E+1  1  , 

T 1 =0 . 0C0 1 524 ,  T2=0. 0001016,  SEND 


*  e 


SEND 


C 


A 
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ft  1) 


CONCLUSION 


The  successful  completion  of  this  Phase  I  study  has  presented  an 
advanced  analytical  tool  for  predicting  the  performance  of  foil  bearings 
and  solid-walled  bearings  as  well,  operating  at  speeds  wherein  the 
Reynolds  number  is  no  longer  small  and  the  effects  of  the  fluid  inertia, 
'  the  heat  transfer  and  the  turbulent  flow  become  significant.  The  tool 
utilizes  the  three  dimensional  Navier-Stokes  approach  for  the  fluid 
film  thermohydrodynamics  and  the  thermoelasticity  treatment  for  the  foil 
deformat i on . 

The  Navier-Stokes  approach  is  much  more  complicated  than  the  Reynolds 
equation  approch,  conventional  or  modified  for  the  inclusion  of  inertia 
or  turbulence.  Nevertheless,  the  N-S  approach  has  the  capability  of 
investigating  the  local  behaviors  of  the  lubricant  and  the  foi'. 
Especially,  the  prediction  of  the  temperature  variation  may  event-ally 
help  explain  the  hot  spots  or  burned  spots  which  are  often  observed  in 
a  fail ed  bearing . 

Naturally,  to  shed  light  on  the  success  or  the  failure  of  a  bearing 
design,  the  steady-state  solution  even  as  complex  as  the  Navier-Stokes 
approach,  -is  just  the  first  step.  Further  studies  should  include  the 
transient  case,  the  combined  bearing-rotor  dynamics  and  possibly  an 
expert  system  for  optimizing  the  design  process. 

The  accomplishment  of  the  present  study  has  established  a  foundation 
for  pursuing  the  just  mentioned  further  steps.  The  formulations,  the 
coordinate  transformations  and  the  numerical  solution  algorithm  car.  be 
extended  tc  deal  with  the  transient  case  and  to  incorporate  with  the 
rotor  dynamic  analysis. 

Two  computer  programs,  3DNSBRGJ  (for  the  journal  bearing)  and  2GNS3PG7 
)  (for  the  thrust  bearing)  are  developed  through  this  Phase  I  study.  These 
rcgrams  are  capable  for  predicting  the  bearing  performance  characteris- 
ics  and  calculating  the  velocity  components,  the  temperature  distri¬ 
bution,  the  foil  deflection,  the  turbulence  energy  distribution  and  its 
dissipation  rate. 
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